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We present an experimentally realizable stabilized charge pumping scheme in a linear array of 
Cooper-pair boxes. The system design intrinsically protects the pumping mechanism from severe 
errors, especially current reversal and spontaneous charge excitation. The quantum Zeno effect is 
implemented to further diminish pumping errors. The characteristics of this scheme are considered 
from the perspective of improving the current standard. Such an improvement bears relevence to 
the closure of the so-called measurement triangle (see D. Averin [Nature 434, 285 (2005)]). 
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Adiabatic variation of parameters in a quantum sys- 
tem results in evolution in which population of each of 
the energy levels is conserved (adiabatic following). It is 
well known that, during such evolution, the energy levels 
acquire two types of phases: time dependent dynamical 
phases and time independent geometrical phases deter- 
mined only by the geometry of the path traversed in pa- 
rameter space. The property of adiabatic following can 
be also used to induce transport through a system, viz. 
a slowly transported potential well carries its bounded 
particles. 

For closed paths, the initial and final states (usually) 
differ only by a phase factor (HI, but the system evo- 
lution during the cyclic variation can still be character- 
ized by non- vanishing transport properties. This phe- 
nomenon, known as parametric pumping, is observable 
in mezoscopic devices, where during each pumping cy- 
cle a definite charge is transported through the system 
P, Q, 0, 0, [H, @l ■ Quantum pumps are hoped to trans- 
form an AC signal of frequency / into DC current given 
by the relation / = fQ, where Q = ne for normal de- 
vices, and Q = 2ne for superconducting pumps (n is 
the total charge transported per cycle). If operated with 
perfect accuracy, quantum pumps could be utilized to es- 
tablish a standard of current. This is needed to close the 
measurement triangle relating voltage, current and fre- 
quency by fundamental constants Q- However, nonadia- 
batic corrections, uncontrolled tunneling, and sensitivity 
of the devices to parametric fluctuations smear the out- 
put signal. An intristic stabilizing property, like quan- 
tization of magnetic flux, crucial for the voltage stan- 
dard, is absent in the parametric charge pumps. Thus 
the only way to stabilize the performance is to optimize 
the design of such devices. E.g. in Ref. [H enhanced con- 
trol over the tunclling amplitudes in the so-called Cooper 
pair "sluice" resulted in improved characteristics of the 
pumping scheme. 

In this Letter, we combine advantages of the control- 
lable Josephson couplings used in the "sluice" and the 
array-like design of earlier proposals [3, 0] . Our scheme 
is analyzed in the framework of adiabatic passage in 



an effective two-level system. Electrostatic coupling of 
the Josephson junctions forming the arrays results in in- 
creased separation of the operational subspace from the 
excited states. The array design also suppresses the most 
severe process affecting the pumping accuracy - current 
reversal. This suppression increases exponentially with 
the number of junctions in the system. Since after each 
pumping sequence the system is to be found in a known 
charge state, we also make use of the quantum Zeno effect 
to diminish the population of unwanted states. 

The core mechanism of the adiabatic passage can be 
understood using the simplest nontrivial - two-level sys- 
tem. Suppose that we start at degeneracy (zero field in 
the spin-1/2 language) and spin pointing "up" in the z 
direction. Then the field is increased along the spin di- 
rection and traverses a path in the x — z plane which 
approaches the degeneracy again from the negative z di- 
rection. For adiabatic pulses, the system remains in the 
ground state and follows the field. However, at the end 
the spin points in the negative z direction. This simple 
example shows that although the system is in a nondc- 
generate level during the passage and the traversed path 
in the parameter space is closed the resulting transforma- 
tion is highly nontrivial as the spin direction is reversed. 
This is due to the level crossing occurring at the begin- 
ning and end of the path. For zero field the states are 
degenerate, and this freedom of defining the orthogonal 
basis results in discontinuity of the energy eigenstates. 

A Cooper-pair box with tunable Josephson cou- 
pling 0] [see Fig. QJa)], described by the Hamiltonian 
H = -l/2B z a z - l/2B x a x , where B x = J($) = 
2 Jo cos(7r<I>/(f>o), and B z = 4Ec(l — 2n g ), mimics a spin- 
1/2 system. Here Jo is the Josephson coupling of the 
junctions in the SQUID (Superconducting Quantum In- 
terference Device) coupling the box to the reservoir, <f> 
and <I>o the external applied flux through the SQUID 



and flux quantum respectively, Ec 



l /(2C g + 2Cj) 



is the single-electron charging energy of the box and 
n g = C g Vg/2e is the dimensionless gate charge. The two- 
level approximation holds if the system is operated in the 
charge regime (Jo <C Ec) and n g is close to 1/2. Then 
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FIG. 1: (a) Cooper pair box[j|. Electrostatic potential of the 
superconducting island (black node) is controlled by the gate 
voltage V g . The Josephson tunelling amplitude J is tuned by 
external flux "1?. (b) A Cooper-pair box with two SQUIDS. 



the only relevant states are characterized by zero or one 
Cooper pair in the box, corresponding to spin parallel 
and antiparallel to the z direction respectively. Perform- 
ing the procedure described above on the system, we can 
thus transport a Cooper pair to or from a reservoir into 
an empty box, depending on the sequence of pulses. 

To induce transport through the Josephson-junction 
system we use two SQUIDs which now serve as terminals 
[see FigQJb)]. Using one of the terminals to transport 
a Cooper pair into the box (stage one) and another to 
transport the same pair to the reservoir (stage two) we 
generate current / = 2e/. Exemplary pulses of a full 
cycle are depicted in Fig. [2j No transport takes place 
during parts I and IV with the system in a definite (zero) 
charge state, so we can limit the cycle to the sequences 
II and III only. This simplification minimizes the cycle 
duration and eliminates degeneracies from the procedure 
(this system is identical to the Cooper pair "sluice" [f|). 

Let us comment briefly on the possible sources of im- 
perfections in the pump performance. The nonadiabatic 
corrections leave the system in an unknown superposition 
of charge states after the full cycles (instead of a definite 
charge state). These corrections whilst usually small, can 
accumulate over time. Fluctuations of external param- 
eters affect the scheme similarly, though fluctuations in 
external flux have a much more severe effect than those 
in the gate voltage, especially when the SQUIDs are not 
ideally symmetric. The terminals cannot then be com- 
pletely closed, due to residual Josephson coupling, and 
the fluctuations induce transitions between charge states. 
The error related to this feature is the hardest to elim- 
inate - since there is nonzero tunneling to the reservoir 
via both terminals all the time. Cooper pairs can acci- 
dentally be transported in the wrong direction, leaving 
no trace of the error in cases when the the final charge 
state is identical to the expected state. In the follow- 
ing we show that the array design combined with Zeno 
projections strongly suppresses these pumping errors. 

The system we consider is shown in Fig. [3ja). It is 
composed of TV Cooper-pair boxes. Nearest neighbours 
are coupled with dc-SQUIDs, leftmost and rightmost is- 
lands are coupled to the reservoirs via other SQUIDs 
(terminals). If the number of Cooper pairs on the is- 




FIG. 2: A pulse sequence applied to the simple pump that 
induces flow of one Cooper pair through the system. The 
exact shape of, and area under the pulses do not influence 
the process as long as the adiabaticity condition holds. 
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, n^v}, the electrostatic energy of the 
= Ei E c(i)(ni-n g(i) ) 2 +E% +1 (ni- 



n 9 (i))(n i+ i - %( i+ i)). Here E c ^ are the charging ener- 
gies of the islands and E^ +1 is the energy of electrostatic 
coupling between two neighboring islands, which is finite 
for finite capacitance of the coupling SQUID. If the di- 
mensionless gate charges, n g u-\ = C g ^V g ^/2e are close 
to 1/2, the lowest energy states are characterized by ei- 
ther zero or one Cooper pair on each island. With this 
assumption we can again reduce the Hilbert space and 
map our system to a finite anisotropic Heisenbcrg spin- 
1/2 chain in an external magnetic field. The Hamiltonian 

1 1 1 N 

H = --Biol ±B»a% - ± ]T Bl4 (1) 

i=l 

N-1 

?:=i 

is characterized by constant electrostatic coupling ampli- 
tudes and tunable parameters - B^' N - the Joseph- 
son coupling of the leftmost and rightmost SQUIDs, B\ - 
electrostatic potential of the islands and Ji,i+\ - Joseph- 
son coupling between neighbouring islands. In general 
the electrostatic coupling has a finite range determined 
by the screening length A = \J Cj / C g (Tcjj . As discussed 
below, we consider the regime in which the electrostatic 
coupling is much weaker than the on-site electrostatic 
energy. This implies that A 1 and we can thus limit 
ourselves to nearest neighbour electrostatic coupling only. 

Because of the electrostatic interaction, for n l g close to 
1/2 the most favorable charge states are antiferromag- 
netic (|010101 . . .) and |101010 . . .)). The pumping pro- 
cedure is a simple generalization of the scheme described 
for one island and two terminals. We start with one of 
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FIG. 3: (a) The stabilized Cooper-pair pump. It is operated 
similarly to the single-island system - identical pulse sequence 
is applied to every second junction (island) (b). Single elec- 
tron transistors (SET) serve as the dephasing elements. They 
are biased with short voltage pulses every time the pump 
should be in a definite charge state [lowest plot in (b)]. 



the antiferromagnetic states and transfer the charge of 
every island to the right by two sites to achieve pump- 
ing. This can be implemented by applying identical pulse 
sequence to every second island (and SQUID), i.e. B\ = 
J2.3 = J4.5 — ■ ■ ■ = Jeo, Ji,2 = J3A = ■ ■ ■ = B x = J oe , 
Bl = Bf = ... = B° and B 2 = Bf = ... = Bf [see 

Z Z Z ' Z Z Z L 

Fig.OJb)] (o, e denote odd and even islands respectively). 
Here we focus on systems with odd number of islands - 
then the number of coupling SQUIDs is even, and exactly 
half of them are active during each step of the procedure. 
This symmetry simplifies our analysis but is not crucial. 

The procedure, even for the single-island realization, 
takes place in nondegenerate ground state. Due to 
the avoided degeneracy, spontaneous transitions to un- 
wanted levels cost energy and are thus suppressed. The 
chain design additionally provides the following advan- 
tage: on the one hand an energy gap between the anti- 
ferromagnetic states and the remaining subspace AE > 
min{Aj.j+i}, while on the other hand accidental transi- 
tions between the antiferromagnetic states need simulta- 
neous transport through (N+l)/2 junctions - another ef- 
fect which diminishes exponentially with TV (see Fig. 2th) 
and explanation below). The most evident advantage 
of this approach is, however, elimination of transport in 
the wrong direction. Let us analyze the mechanism rul- 



ing this process. Suppose that we start with the state 
1 101010 . . .) . During the first part of the cycle we want to 
shift the charge configuration by one site to the right, i.e. 
to arrive at |010101 . . .). As mentioned before, the resid- 
ual Josephson coupling of closed SQUIDs leads to two 
ways of arriving at the same state. One could anticipate 
interference effects here, but since the residual coupling 
is a fluctuating quantity (controlled by fluctuating mag- 
netic field) we treat the process as incoherent. A simple 
analysis is based on the assumption that a Cooper pair 
is transported through an open (correct) SQUID with 
probability p, and through closed (wrong) SQUID with 
the probability q — 1 — p. The measure of this pump- 
ing error for each island might be the ratio q/p. For 
a chain of N junctions the probability of correct shift 
becomes oc p N ^ 2 , and of incorrect shift oc q N l 2 . Hence 
the pumping error (q/p) N ^ 2 decreases exponentially with 
the number of junctions in the chain. More rigorous 
analysis requires calculation of the energy difference be- 
tween the two lowest energy levels for B\ = 0. Indeed, 
in this case the lowest energy eigenstates are symmetric 
and antisymmetric superpositions of the antiferromag- 
netic states (with some addition of other charge states), 
and the energy difference, being the rate of phase pro- 
gression, in the charge basis is equivalent to the rate of 
transition. Fig. [4f a) shows the results for various val- 
ues of the maximum Josephson coupling of the SQUIDs. 
T is the rate calculated for the desired transitions only 
(assuming no residual coupling), t is calculated for the 
wrong path. Then, the time of single operations should 
be much larger than 1/T for the adiabatic condition to be 
satisfied. The numerical calculation is performed for a re- 
alistically achievable residual coupling equal to 1% of the 
maximum Josephson coupling of the SQUIDs. The values 
of the electrostatic coupling amplitudes, and maximum 
Josephson couplings are for simplicity assumed identi- 
cal for each junction (this is, again, not crucial and our 
observations are valid also for different settings). The 
pumping error vs. N is plotted in Fig.[4]^a). The depen- 
dence is clearly exponential, and very weakly depends 
on the value of the maximum Josephson coupling. How- 
ever, the value of T itself should not decrease too strongly 
with increasing N as it determines the adiabatic rates. 
From Fig.[4|b) we can see that the ratio J max /A should 
be as close to 1 as possible to avoid the effects of ex- 
ponential suppression. For this case, T/A « 1, where 
A can be chosen large thus yielding a very comfortable 
adiabatic condition. Note that even for J max /A sw 1 
we can still work in the reduced Hilbert space, as long 
as the charging energies of individual islands are much 
larger than J max (this means that the system should be 
operated in the charging regime with weak electrostatic 
coupling). The exponential suppression of the energy gap 
between the lowest states implies that unwanted transi- 
tions between the antiferromagnetic states also diminish 
very rapidly with increasing N (note that this suppres- 
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FIG. 4: (a) Pumping error t/T as a function of N. For log- 
arithmic vertical axis the dependence is linear, implying ex- 
ponential suppression of the error with increasing TV. (b) 
Transition rates between the charge states as a function of N. 



sion for very small, erroneous, values of Josephson cou- 
pling is much stronger than for moderate values, expected 
for open terminals). This of course does not mean that 
the optimal setup should be represented by as long an 
array as possible. Most probably there is some moder- 
ate number of junctions for which the pumping errors 
are strongly suppressed, and the decoherence effects typ- 
ical for macroscopic system are yet not severe. However, 
examination of the long-chain limit goes beyond our anal- 
ysis and would be appropriate for a separate publication. 

Although during the passage the state of the system is 
a superposition of numerous charge configurations, when 
all Josephson couplings are switched off (or rather in- 
tended to be) the pump should be in one of the an- 
tifcrromagnctic states. All other components (nonadi- 
abatic corrections, residual Josephson coupling effects) 
contribute to pumping errors. The optimal solution 
would be to project the state onto a charge state by per- 
forming charge measurement, fndced, if the measured 
state is very close to an eigenstate of the measured ob- 
servable, the measurement takes the state closer to the 
eigenstate. If frequently repeated, the measurement can 
diminish the probability of transition to other eigenstates 
arbitrarily. This phenomenon, known as quantum Zeno 
effect , can be utilized also in our system. The charge 
measurement can be performed by a single-electron tran- 



sistor (SET) coupled capacitively to the superconducting 
island [l2[ [see Fig. [3ja)]. This measurement is charac- 
terized by few time scales, the dephasing time r^, during 
which the off-diagonal density matrix elements (in the 
basis of the measured observable) vanish, the measure- 
ment time, T measi after which the information about the 
charge can be read out, and the mixing time, after which 
the back action of the SET on the island destroys the 
information about the state. The parameters can be se- 
lected in such way that r m „ > r meas > tq. Since the 
actual information about the state is not important, and 
only the dephasing mechanism modifies the state, the 
quantum Zeno projection can be realized within a rela- 
tively short time. Moreover, even if the time available is 
too short for the off-diagonal entries of the density ma- 
trix to vanish, partial dephasing still brings the system 
closer to the desired charge state. The system can be 
thus dephased each time the system is supposed to be in 
a definite charge state by switching on the bias voltage 
of the SETs for a very short time interval. 

In summary, the presented scheme provides an accu- 
rate pumping procedure. The array design exponen- 
tially suppresses the most severe, untracable error of re- 
versed current flow and diminishes errors due to unde- 
sired charge configurations. These errors can be further 
minimized by repeated dephasing with single-electron 
transistors. 
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